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Abstract 

Corrections of 0{al) to the decay of the top quark into a W boson and a bottom 
quark are calculated. The method is based on an expansion of the top quark propa- 
gator for small external momentum, q, as compared to the top quark mass, Mj. The 
physical point q^ = Mf is reached through Pade approximations. The described 
method allows to take effects induced by a finite W boson mass into account. The 
numerical relevance of the result is discussed. Important cross-checks against recent 
results for the decay rate 6 — > u/P and the two-loop QED corrections to /i-decay are 
performed. 
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1 Introduction 

The top quark is the so far heaviest observed particle of the Standard Model of elementary 
particle physics. Its total width, Tt, is to a good approximation proportional to the third 
power of its mass and is much larger than Aqcd, the typical scale of non-perturbative 
effects in QCD. Therefore it is possible to treat the top quark almost as a free particle 
and to apply perturbative methods to evaluate the quantum corrections to its decay 
process |]I]]. 

In the minimal Standard Model the dominant decay mode of the top quark is the one into 
a bottom quark and a W boson. It is important to predict the corresponding decay width 
accurately in order to be sensitive to exotic processes, which may occur in supersymmetric 
models, for example. 

The first order QCD corrections have been evaluated in analytical form some time ago 
and amount to approximately —10%. The electroweak corrections are small [^ and turn 
out to be ^ 2% for a Higgs mass around 100 GeV (see, e.g., 0). 

The expected precision for measurements of Tt by a future e~^e~ machine like the Next- 
Linear-Collider (NLC) is of the same order of magnitude as the corrections of 0{as) ^■ 
This makes it desirable to control also the next-to-leading corrections induced by the 
strong interaction. 

In fact, the QCD corrections of 0{al) have already been considered in [^]. This calculation 
was based on an expansion of the vertex diagrams in the quantities 1 — M^/M^ and 
1 — SM^/M"^, respectively. Although this expansion parameter is not small at all, the 
approach led to reliable results after including many terms into the analysis, choosing 
proper variables, and carefully investigating potentially large contributions. 

The aim of this paper is, on the one hand, to have an independent check of the results 
of , using a rather complementary method. On the other hand, our approach will allow 
us to additionally account for a finite Myy boson mass. 

The method presented in this paper is as follows. In contrast to |^ we compute propaga- 
tor-type diagrams contributing to the top quark selfenergy with external momentum q in 
terms of an expansion around q^ /Mf = 0. Some sample diagrams are pictured in Fig. [^. 
The imaginary part combined with the wave function renormalization of the top quark 
and evaluated at the physical point g^ = M^ directly leads to the decay rate. It arises 
from cuts where the W boson, the bottom quark and, at higher orders, also gluons and 
other light quarks are involved. The limit q^ —>■ M^ is taken after performing a Fade 
approximation. The results for M^ = will be shown to be in perfect agreement with 
the ones of which justifies both the method of ||^ and the one of the present paper. 
The subleading terms in M^ turn out to be numerically small. 

The calculation once again demonstrates the power of expansion techniques and their 
computer implementations in multi-loop calculations. The analyticity properties of the 
approximated function guarantee reasonable convergence to the exact result, especially if 



w w 

b\ t t Jb 



w w w w 

b\ t t Jb\ t t Jb\ t t Jb 





u/d/s/c/b t 

Figure 1: Sample diagrams for the top quark self energy. 

the obtained series is further subject to advanced methods like Fade approximation. 

The paper is organized as follows: In Section ^ the method is described. In Section ^, the 
results obtained at 0{as) are discussed in more detail. The comparison with the analytical 
result demonstrates the reliability of our method. Section ^ deals with the computation 
of the second order QCD corrections where also the effects of a finite W boson mass 
are taken into account. Important cross-checks with recent results for F(6 -^ ulv) and 



F(/i -^ eUf^Ue) are carried out in Section ^]5 



2 Method of the calculation 

The exact evaluation of the Feynman diagrams contributing to r(t -^ Wb) at order a^ 
is currently not available. However, it is promising to apply the methods of asymptotic 
expansion (see, e.g., and references therein) in the limit 



M^ > g2 > M^ > M^ = . 



For the 0{as) corrections in the limit M\y = 0, a similar approach has already been 
used in [0]. There it was possible to resum the series at g^ = M^ which reproduces the 
analytical result. At 0{a1), however, instead of an explicit resummation we will perform 
a Fade approximation in order to reach the physical point g^ = M^ [§, ^. 

Before going into details, let us introduce the notation. The inverse quark propagator is 
denoted by 



z mMl-S^ -^ l + S"^ . (2) 



Both S5. and S^ are functions of the external momentum q and the bare mass m^ of the 
top quark. In our case they further depend on the bare strong coupling constant a° and 
the W boson mass M^, and are proportional to the Fermi coupling constant, Gp- Sp 



will become finite after renormalizing the parameters m^ = Z^^Mf and a^ = Z^ag, and 
taking into account the wave function renormalization: 

^F ~ '70S F- (3) 

Z^'^, Z2^ and Zg denote the renormalization constants. Z^^ and Z^"^ will be taken in 
the on-shell scheme, whereas Zg is in the MS scheme. Z^*^ is defined by the condition 

Sf^{q) — ^ 1- + terms regular for q^ -^ M'f . (4) 

Mt - i 

In our approach we are actually dealing with two different masses for the top quark in 
intermediate steps: an "external" one (g^) and an "internal" one [M^). Applying the 
optical theorem, the decay rate of the top quark will therefore be written as 



T{t -^ Wb) = {2Mt lm[z ^^^ - 5^^]) 
where 



(5) 



are functions of the variable 
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(6) 



(7) 



All relevant diagrams will be calculated in terms of expansions around z = 0, and the 
limit z ^ 1 will be applied only in the very end. Therefore, while Z^^ can be taken at 
2; = 1, we also need to express Z^*^ in terms of an expansion around z = 0. This is most 
conveniently done by translating condition (HI) into 



Z2 



OS 



l + S°+2^fz'^^S°-S° 



.y^.-y^^^S ^V 



From (^ and (P) it is clear that for our purpose it suffices to know only the pure QCD 
corrections of Z^"^ and Z^^ up to O^a"^). For z = 1 these quantities were computed up 
to this order in [^ and |T^, respectively. 



Note that in a calculation where the quantities are evaluated on-shell, i.e. at g^ = M^, 
infra-red singularities appear in intermediate steps. In contrast, in Eq. (|^) all functions 
on the r.h.s. are defined through the expansion for 2; — > and thus are infra-red safe. 

At this point a comment on the extraction of the values for z = 1 is in order. Actually 
Eqs. (H) and (H) are not unique as it is possible to derive slightly different equations for 
Z2^ and r(t ^ Wb), which differ by relative factors of 2;. In the limit z ^ 1 all of them 



are equivalent. The results we obtained by using two more variants of Eqs. (^) and (P) are 
consistent with the ones which will be discussed below. We decided to use the formulae 
shown above because the 0{as) corrections are recovered with the highest accuracy. 

In order to obtain reliable results it is necessary to evaluate as many terms as possible 
in the expansion parameter z. The exact resummation of the series in z seems to be 
excluded. Instead, we apply a Pade approximation which means that we reexpress the 
resulting polynomial in terms of a rational function: 



^m 



ao + aiz + . . . + a^z ' , , 

\m/n\[z) = — ; ; . (9) 

Its Taylor series is required to coincide with the original polynomial up to the order m + n. 
For later convenience we define the short hand notation [m/n] = [m/n]{l). The stability 
of the Pade approximants upon variation of m and n will indicate the uncertainty of the 
approximation (see below). 

In addition, it may be promising to apply a conformal mapping P] 

4a; 

(1 +UJ)^ 

and to perform the Pade approximation in the variable u. The complex z-plane is mapped 
into the interior of the unit circle in the u;-plane, and the relevant point z = 1 goes to 
uj = 1. This conformal mapping is motivated by the observation that the application of 
a Pade approximation relies heavily on analytic properties. The function on the r.h.s. of 
Eq. @ (without the limit z —>■ 1) will develop a branch cut along the real ^-axis starting 
from z = 1. This branch cut is mapped through Eq. (^) onto the unit circle. Thus by 
applying Eq. ([T^) we enlarge the range of convergence for the terms we got in the limit 
z-^0. 

Since both methods described above appear to be rather natural, any of them will be 
used to derive an estimate on the exact result. For convenience, let us denote the results 
obtained through Pade approximation in the variable z by "2;-Pades", the ones where 
the Pade approximation is performed in uj by "cj-Pades". The central values and the 
estimated uncertainty will be extracted from Pade results [m/n] with m + n not too 
small and \m — n\ < 2. The central value is obtained by averaging the Pade results and 
the uncertainty is given by the maximum deviation from the central value. The error 
estimation is therefore rather conservative. 

Some Pade approximants develop poles inside the unit circle {\z\ < 1 and \uj\ < 1, 
respectively). In general we will discard such results in the following. In some cases, 
however, the pole coincides with a zero of the numerator up to several digits accuracy, 
and these Pade approximations will be taken into account for the estimation of the actual 
results. To be precise: in addition to the Pade results without any poles inside the unit 
circle, we will use the ones where the poles are accompanied by zeros within a circle of 
radius 0.01, and the distance between the pole and the physically relevant point q^ /Mf = 1 
is larger than 0.1. 



Concerning the dependence on the strong gauge parameter ^ in Eq. (^ , it only drops out 
after summing infinitely many terms in the expansion around z = and setting z = 1. 
Since we are only dealing with a limited number of terms, our approximate results will 
still depend on the choice of ^ even after taking z ^ 1. It is clear that with extreme 
values of ^ almost any number could be produced. Thus the question arises which value 
of ^ should be assumed in order to arrive at a reliable prediction for r{t ^ Wb). 

At 0{as) the whole calculation can be performed for arbitrary gauge parameter without 
any difficulties. This allows for a detailed study of the residual ^ dependence. At C(a^) 
only the first few terms could be evaluated for general ^ which does not allow for extensive 
studies. In order to arrive at a reasonable number of terms in the expansion around z = 
it is necessary to set ^ to some definite value from the very beginning. The behaviour of 
the diagrams at 0{as) will serve as an indication for the optimal choice of this value in 
the analysis at 0{a1). 

Concerning the electroweak gauge parameter ^w, all results that will be quoted in the 
following have been obtained in unitary gauge, where the W propagator is given by 

Nevertheless, the leading terms in Mw have also been computed in an arbitrary covariant 
gauge. They are obtained by replacing the W boson by a Goldstone boson with the 
propagator simply given by 

The independence of C,w is then manifest already at this point. 



3 First order QCD corrections 

In this section we will investigate the 0{as) corrections and compare the exact result 
to the approximation obtained by the method described above. 

It is convenient to decompose the decay rate of the top quark into a W boson and a 
bottom quark in the following way: 



r(t ^ bW) 



a;, „ .n^ fa^^"^ 



TT V TT 



(13) 



where Tq = GFMf\Vtb\'^ /{^ti^/2), A^ = 1 - 3M^/M/ + 2M^/Mf, Cp = 4/3, and Vtb is 
the CKM matrix element for t ^ b transitions. The running coupling a^ is defined with 
six active fiavours. 



The one-loop correction is known in analytical form since quite some time 
in terms of Myy/Mt it reads: 
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(14) 



The approximation Mw = induces an error of roughly 22%. This reduces to approxi- 
mately 4% if the quadratic mass corrections are included and is completely negligible if 
all the terms displayed in Eq. (|1^) are taken into account. 

For clarity, let us apply our method to these lowest order terms and see how the results 
compare to A'^'^'^ and A^^^ above. While A^^^ is reproduced exactly, the imaginary part of 
the small-momentum expansion for the two-loop {0{as)) diagrams (an example is shown 
in Fig. 0) reads: 
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(15) 



with Imi = In M^/M^. (The coefficient of z^ will be called the "n*'' moment" in the 
following.) Note that the factors 1/z and l/z"^ in front of the quadratic and quartic terms 
in Mw are irrelevant for the subsequent Fade procedure. 



The procedure described above is applied to each coefficient of M'yy/M^ separately. As 
already noted, A'^^^^ still depends on the QCD gauge parameter, ^, appearing in the gluon 
propagator i{—g^^ + S,q^q^ /q^)/{q^ + ie). Thus also the Pade approximations [m/n] will 
show a dependence on ^. It is clear that for large absolute values of ^ the quantities [m/n] 
get dominated by them and any predictive power is lost. In Table |l| several z-Pades are 
evaluated for the leading order coefficient {Mw = 0). The gauge parameter is varied from 
^ = —2 to ^ = +2. Pade results which develop poles for |2;| < 1 are in general represented 
by a dash. However, if an approximate cancellation with a zero from the numerator 
takes place (see the discussion above), they are marked by a star {-k). The exact result 
is reproduced with a fairly high accuracy for almost all values of ^ under consideration. 
Nevertheless, the value for ^ = is closest to the exact result, and the variation of the 
Pade approximants appears to be very small for this particular choice of ^. Based on this 
observation, we decide to perform the three- loop analysis by setting ^ = 0. This has the 
additional advantage that the evaluation of the Feynman diagrams is much simpler than 
for non-zero values of ^. 



input 


P.A. 


e = -2 


e = -l 


e = -l/2 


e = o 


e = l/2 


e = i 


e = 2 


6 
6 


[3/2] 
[2/3] 


-2.111 
-2.112 


-2.052 


-2.058 
-2.058 


-2.023 
-2.023 


-1.990 
-1.990 


-1.957 
-1.963 


-1.893 

-1.887 


7 
7 
7 


[4/2] 
[3/3] 
[2/4] 


-2.121 
-2.120 
-2.126 


— 


-2.058 
-2.058 
-2.058 


-2.027 

W - 2.025 

-2.027 


-1.999 
-2.008 
-2.000 


-1.972 
-1.980 
-1.979 


-1.919 
-1.928 
-1.912 


8 

8 


[4/3] 
[3/4] 


-2.117 
-2.117 


-1.694 


-2.058 


-2.033 
-2.033 


-2.011 
-2.011 


-1.990 
-1.993 


-1.949 
-1.948 


9 
9 
9 


[5/3] 
[4/4] 
[3/5] 


-2.112 
-2.105 
-2.114 


-2.063 
-2.061 
-2.064 


(*) - 2.059 
W - 2.059 
W - 2.059 


-2.034 
-2.034 
-2.034 


-2.016 
-2.040 
-2.016 


-1.998 
-2.001 
-2.001 


-1.963 
-1.967 
-1.961 


exact: 


-2.040 



Table 1: z-Pades for A^^) 



M„ 



for different values of ^. The first column 



indicates the number of terms in z from (^) that were used as input. 



Let us now discuss the sub-leading terms in M^. Table |^ lists several Pade approximations 
for the coefficients of (M^/M^)" [n = 0, 1, 2) in the case of A^^\ ^ has been set to zero, 
as it is suggested by the discussion above. For each coefficient, the Pade approximations 
have been performed in the variable z as well as in the variable u. The z- and tu-Pades are 
indicated by a z and an u, respectively, in the second line of Table |^. This notation will be 
used throughout the paper. The z-Pades for the values of the {M^r / M^)^ -term coincide 
with those for ^ = of Table |I], of course. Concerning the power-suppressed terms, 
again the higher order Pade approximants agree with the exact results to an impressive 
accuracy. The logarithm of Eq. (|I4D is reproduced exactly after taking into account the 
ffist two terms in the expansion (|I^) and setting z = 1. 

Also for the power-suppressed terms an analysis concerning the ^ dependence has been 
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to 
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LO 


6 
6 


[3/2] 
[2/3] 


-2.023 
-2.023 


-2.022 
-2.051 


1.506 


1.508 
1.548 


3.875 

3.878 


I 


7 
7 
7 


[4/2] 
[3/3] 
[2/4] 


-2.027 

(*) - 2.025 

-2.027 


(*) 


- 2.009 
-2.035 
-2.040 


Wl.507 


(*U.518 
1.501 
1.507 


3.871 
3.872 
3.874 


— 


8 
8 


[4/3] 
[3/4] 


-2.033 
-2.033 


-2.035 
-2.035 


1.502 


1.502 
1.502 


3.870 
3.870 


I 


9 
9 
9 


[5/3] 
[4/4] 
[3/5] 


-2.034 
-2.034 
-2.034 


(*) 


-2.039 
- 2.035 
-2.037 


1.502 
1.502 
1.502 


1.500 

Wl.501 

1.499 


3.870 
3.870 
3.870 


— 


exact: 


-2.040 


1.500 


3.870 



Table 2: Pade approximations for the power-suppressed terms of y4^^\ com- 
puted for .^ = 0. 



performed. While the quadratic terms in Ms^ do not depend on ^ at all (see Eq. (|T5|) ), 
the results for the terms of order M^jM^ become significantly worse once ^ is different 
from zero. This also supports the choice of Feynman gauge at three loops. 

Taking only those results of Table |^ into account where eight or more input terms enter 
we may finally deduce our approximation for the 0[a^ corrections to the decay rate (we 
adopt the notation -2.035(4) = -2.035 ± 0.004, etc.): 



A(^) = -2.035(4) 
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The agreement of Eq. (0) with the exact results quoted in the last line of Table [ 
supports the underlying method. 



I obviously 



In the three-loop case we could obtain the small-momentum expansion up to z^ for My^ = 
and up to z^ for the coefficients of M^/M^^ and the M^r/Mf. Therefore, the final number 
at three-loop level will be based on Fades built out of seven and eight moments for the 
leading term, and six and seven moments for the sub- leading terms in My/. At 0{as) 
this reduced number of input terms changes the result from the one in (fB) to 
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which is a bit worse than ([T6|), but still sufficiently accurate. This suggests that the 
number of available moments at order a^ should be sufficient to arrive at a reasonable 
estimate. 



4 Second order QCD corrections 

Let us use the experience gained in the previous section to obtain predictions for T{t 
Wb) at order a^. 



4.1 General remarks 

It is convenient to decompose the decay rate according to the colour structure: 

^(2) = C^pA^J^+CACpA^NA + CFTniAi^^+CFTAP, (18) 

where in QCD the colour factors are given by Cp = 4/3, Ca = 3, and T = 1/2. rii is the 

(2) 

number of massless quark flavours and will be set to n; = 5 in the end. A\ corresponds 

(2) 

to the abelian part already present in QED, A)^]^ represents the non-abelian contribution, 

(2) (2) 

and A] and Ap denote the corrections involving a second fermion loop with massless 
and massive quarks, respectively. In Fig. |I| a representative diagram for each function is 
pictured. The expansion in terms of M^/M^ of the individual contributions to A^"^^ will 
be written as 

4(2) _ 4(2), , M^ .(2), M^ (2) 

with i = A, NA, I and F. Note that Ai \mw=o is known analytically [|l^ and serves as a 
welcome check for our method. 



a 



2A 



Whereas at 0{as) the 't Hooft mass /i^ drops out (see Eq. ([T^)), it does appear at 0{ 
We adopted the convention /x^ = M| throughout the paper. 

There are 60 three- loop diagrams that contribute to r(t -^ Wb). The practical com- 
putation is done with the help of the package GEFICOM [|13[. It uses QGRAF ||T4l for the 
generation of the diagrams and EXP [|15| for the application of the hard mass procedure. 
For more details we refer to a recent review concerned with the automatic computation of 
Feynman diagrams . The application of the methods of asymptotic expansion according 
to Eq. (|I|) reduces the practical computation either to massless three-loop propagator-type 
diagrams or to products of one- and two-loop integrals. In the latter case either vacuum 
graphs or again massless two-point functions appear. The integrals have been performed 
with the help of the packages MINCER []T6| and MATAD [ll^ based on the symbolic manip- 



ulation language FORM |T^. The results for S^ and Sy defined in Eqs. (^ and (H) at 
three loops are quite lengthy and therefore not listed here. They can be obtained from 
the authors upon request. 

4.2 The limit Mw = 

This section is concerned with the second order QCD corrections where the mass of the 
W boson is neglected. In this limit a comparison with |^ can be performed. The Mw- 



suppressed corrections will be discussed in the subsequent section. As already noted, 
'^A,expi^) S'^d ANA,expi^) Were computed up to z'^ in the case of a massless W boson in 
Feynman gauge which corresponds to eight input terms for the Fade approximations. 

The fermionic pieces, AlJ^^ and Ap^^^^, do not depend on the QCD gauge parameter and 

(2) (2) 

are of simpler structure than A)^'^^p and ^TVA.exp- For the light-fermion contribution we 
thus could evaluate nine terms in the expansion around z = for the leading term in 

In Table ^ the results are displayed. The z-Fades for Ap converge very quickly whereas 
most of the tu-Fades develop poles for |ct;| < 1. As a result we infer 



A 



(2) 



Mw=0 



-0.06360(1) 



(20) 



l(2) 



which coincides with the one quoted in 0. Note, however, that the magnitude of Ap is 
rather small. 
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Af 


Af 




input 


P.A. 


z 


UJ 


z 


UJ 


z 


UJ 


z 


UJ 


6 
6 


3/2 
2/3 


2.600 
2.600 


2.889 
2.994 


-0.06359 
-0.06359 


-0.06464 


2.971 
2.972 





-7.637 
-7.639 


-8.238 
-8.400 


7 
7 
7 


4/2 
[3/3] 
[2/4] 


2.633 

W2.592 
2.636 


2.920 
2.926 
2.956 


-0.06360 
-0.06360 
-0.06359 


— 


3.041 
3.140 
3.046 


3.818 


-7.710 
-7.781 
-7.719 


-8.212 
-8.218 
-8.306 


8 
8 


[4/3] 
3/4 


2.695 
2.696 


W2.907 
2.832 


-0.06360 
-0.06360 


— 


3.146 
3.146 





-7.820 
-7.820 


W - 8.232 
-8.254 


9 
9 
9 


5/3 
4/4 
[3/5] 


2.708 
2.707 
2.708 


2.881 
2.892 
2.902 




exact: 


2.859 





Table 3: Fade results at 0{a1) for M\y = 0. 



A] behaves similar to the 0{as) corrections. As expected, the more terms of the expan- 
sion in z are included, the better agreement with the exact result is observed. Furthermore, 
both z-Fades and w-Fades lead to compatible numerical values from which, including the 
seventh and eighth moment ( "input 8" and "input 9" in Table |^) , the following result is 
deduced: 



A\ 



(2) 



2.8fr 



Mw=0 



(21) 



The error is around 4% and thus roughly as large as the one in |^ , where the result reads 
2.85(7). Using only the sixth and seventh moment one ends up with 



A 



(2) 



Mw=0 



2.8(2) 



(22) 



10 



where the error is 7%. The result of Eqs. (^) and (22) can also be compared with the 



exact number [121 which reads 2.859 



Let us now turn to the abelian and non-abelian parts. Like in the one-loop case, the 
expansion in z is gauge dependent for these contributions. Motivated by the observations 
of Section |^, the analysis will be performed by setting ^ = from the very beginning. For 
non-zero values of ^ some lower-order Fade approximants will be presented at the end of 
this section. 

As compared to the fermionic contributions, the spread among the 2;-Fades and tu-Fades 
is significantly larger. Moreover, the numbers for the two approaches are less compat- 
ible with each other. Nevertheless, following the previously introduced strategy for the 
extraction of the central value and the error, we obtain 



A 



(2) 



A 



(2) 
NA 



Afw=0 



3.2(6), 
-8.0(3) . 



(23) 
(24) 



The (fairly conservative) errors are larger than the ones of the results in |^, which read 
3.5(2) and —8.10(17), respectively. The numbers, however, are consistent. 

At this point we have confirmed the results of with a completely independent method. 
Our calculation can therefore serve as an important cross check. 

Due to the complexity of the intermediate expressions it was impossible to evaluate eight 
terms in the expansion for small z using a general gauge parameter. In this case, we 
managed to compute the expansions for Jvj( and A^^\ only up to terms of (9(2^). In 
Tables ^ and ^ the gauge parameter is varied between ^ = —2 and ^ = +2. For A\ the 
error bars of Eq. (p3|) are conservative enough to cover even most of the values of Table |^. 

(2) 

The stability of A}^^ against variations of C, is also satisfactory, although many of the 



values of Table |^ are not compatible with 
number of input terms in Table B. 



. However, this could be traced to the low 



input 


P.A. 


e=-2 


e = -l 


e = -i/2 


? = o 


e = l/2 


e = l 


e = 2 


3 

3 
3 


[2/0] 

[1/1] 
[0/2] 


2.790 
2.780 

2.787 


2.425 
2.432 


2.305 
2.685 
2.329 


2.226 
2.508 
2.280 


2.189 
2.518 
2.293 


2.194 
2.632 
2.375 


2.329 
3.254 
2.789 


4 
4 
4 
4 


[3/0] 
[2/1] 
[1/2] 
[0/3] 


2.849 

W2.624 

W2.705 

2.847 


2.501 
3.268 
3.320 
2.517 


2.397 
2.810 
2.810 
2.433 


2.341 
2.713 
2.721 
2.411 


2.332 

2.725 
2.739 
2.455 


2.369 
2.816 
2.827 
2.572 


2.586 
3.201 
3.202 
3.047 



(2) 

Table 4: ^ dependence of A\ \mw= 
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input 


P.A. 


C = -2 


C = -i 


C=-l/2 


e = o 


e = i/2 


C = l 


e = 2 


3 


[2/01 


-7.894 


-6.953 


-6.571 


-6.248 


-5.984 


-5.779 


-5.547 


3 


[1/11 


-8.202 


-7.455 


-7.176 


-6.954 


-6.783 


-6.655 


-6.515 


3 


[0/2] 


-8.078 


-7.243 


-6.924 


-6.669 


-6.476 


-6.340 


-6.219 


4 


[3/01 


-8.117 


-7.262 


-6.915 


-6.623 


-6.385 


-6.202 


-5.997 


4 


[2/11 


-8.335 


-7.670 


-7.417 


-7.212 


-7.051 


-6.930 


-6.802 


4 


[1/21 


-8.346 


-7.690 


-7.440 


-7.236 


-7.075 


-6.956 


-6.834 


4 


[0/3] 


-8.241 


-7.493 


-7.209 


-6.984 


-6.813 


-6.691 


-6.580 



Table 5: 



^ dependence of A]^'j^\mw= 



4.3 Subleading terms in Mw 

Let us now turn to the power-suppressed terms of order (M^/M^^)" (n > 0). Both for 
the 0{M^r/M^) and 0{M^^/Mf) contribution seven terms in the expansion in z could 
be evaluated. 



4.3.1 Quadratic terms in Mw 





A- 


^P 


4(2) 


.(2) 


input 


P.A. 


z 


LJ 


z 


w 


2, 


w 


Z 


w 


6 
6 


[3/2] 
[2/3] 


W - 1.013 


-1.016 
-0.9459 


0.09766 
0.09766 


0.09769 
0.09769 


(*) - 2.721 
-2.747 


— 


3.358 
3.358 


— 


7 
7 
7 


[4/2] 
[3/3] 
[2/4] 


-0.9737 
-0.9758 
-0.9652 


-0.9773 

-0.9676 

-1.122 


0.09766 
0.09766 
0.09766 


0.09764 
0.09764 


W - 2.677 

W - 2.677 

-2.770 


-2.794 


3.353 
3.354 
3.356 


— 



Table 6: Fade results for the O(a^) coefficients at (^(M^/M^^ 



The resulting Fade values for the quadratic M^y terms are listed in Table ||. The light- 
fermion contribution is very stable and the maximal deviation from the central value 
amounts to roughly 10%. Like in the case Mvk = 0, the Fade approximants for Jvp — 
both in 10 and z — exhibit an impressive convergence. 

For the abelian and non-abelian contribution most of the cj-Fades develop poles inside 

\^\ are very smooth, for Jvj^ there are only 
< 1. However, all the poles of the other 2;- Fades 



\U)\ 



< 1. In addition, while the z-Fades for A^j^ 



two of them without poles within \z 

approximately cancel against zeros in the numerator. 

consistent. 



All relevant numbers are highly 



12 



Therefore, the numbers in Table M lead us to the following results: 



A? 


^l 


= -2.73(6), 


4(2) 


M'w 


= 3.356(3), 


aP 


Ml. 


= -1.0(1), 


aP 


^l 


= 0.09766(3). 



(25) 



4.3.2 Quartic terms in Mw 

In Table the Pade approximations for the quartic terms in Mw are listed. The conformal 
mapping seems to spoil the convergence property here, as all u;-Pades develop poles within 
the unit circle. Moreover, in contrast to the constant and quadratic corrections, for the 
M^ terms it turns out that variations of Eqs. (|) and (|D (see the discussion in Section 0) 
lead to results which lie outside the error interval obtained from the numbers of Table |^. 
So the final numbers for the M^ terms should only be considered as estimates on their 
order of magnitude. They will be presented below, with an artificially increased error of 
about 50%. 





A- 


A2) 






input 


P.A. 


z 


LU 




UJ 


z 


bJ 


z 


OJ 


6 
6 


[3/2] 
[2/3] 


-1.290 
-1.291 


I 


0.1474 
WO.1384 




(*)4.606 





2.364 

2.448 




7 
7 
7 


[4/2] 
[3/3] 
[2/4] 


W - 1.268 

W - 1.278 

-1.305 


— 


0.1479 
0.1477 
0.1489 


— 


(*U.460 
W4.46I 


— 


2.489 
2.459 
2.460 


— 



Table 7: Fade results for the 0{al) coefficients at 0{M^/Mf). 



While the leading two terms in Mw do not contain logarithms of My/-, the coefficients 
of M^ develop linear logarithms of M^r/M^. For A^"^^ the coefficient of this logarithm 
is exactly reproduced by the ffist two terms in the Taylor expansion around z = after 
setting z = 1] the higher order terms in z vanish (a similar behaviour for the logarithmic 
terms has already been observed in [|1^). The phenomenon of a truncated series in z for 

According to the discussion above, this 



12) 
this coefficient also appears for AiJ^p and A 



1(2) 

F,exp' 



strongly suggests that the logarithms of Mw are exactly recovered after setting z = 1. 

While for A} one arrives at 7/4 ■ In M^/M^r for the term under consideration, it sums 

(2) 
up to zero for Ap' . 
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For the abelian and non-abelian parts, the series in z for the coefficients of the My^/Mf ■ 
In M^/Myr term does not seem to be truncated, but the Fade analysis turns out to be 
very stable. It can be found in Table ||. 





.(2) 




input 


P.A. 


z 


U! 


z 


LU 


6 
6 


[3/2] 
[2/3] 


0.7811 
0.7829 




-3.614 
-3.614 





7 
7 
7 


[4/2] 
[3/3] 
[2/4] 


0.7334 
0.7747 
0.7594 


— 


-3.626 
-3.618 
-3.622 






Table 8: Fade results for the coefficient of M^/M^ ■ \n{M^/M^) for aJ^ and 



A 



(2) 
NA- 



Finally, the results for the quartic contributions in M^ read: 



.(2) 
^A 



Am 

^NA 



A\ 



(2) 



A 



(2) 



^^ 



M^ 



AC 



Mtv 



4.5(2.2) + 0.7(1) In 



Ml 



2.4(1.2) -3.62(1) In 
0.15(5). 



Ml 



(26) 



4.4 Results at 0{a^,) 



In this subsection we finally present the numerical corrections for F(t -^ Wb) at 0(a'^). 
Simply using the results derived above according to (|T3p and linearly adding the errors 
one would certainly overestimate the total uncertainty. It is more promising to add up 
the expansions in z for different colour structures and to perform the Fade procedure 
afterwards. The corresponding Fade approximations for the case M^ = are shown in 
Table ^. Table contains the M^iz-suppressed terms. The behaviour is similar to the 
approach where the individual colour structures are treated separately. In the case of 
vanishing W boson mass the Fade results both with and without conformal mapping are 
highly stable. For the M^- and M^- terms, however, many of the cj-Fades develop poles 



inside the unit circle. The result we deduce from Tables and 10 reads: 



^(2) 



-16.7(8) + 5.4(4) 



Ml 
Mf 



+ 



M^ 
Mt 



11.4(5.0) -7.3(1) In 



Ml 



(27) 
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The leading-order result is in very good agreement with the one of which reads 
—16.7(5). Note again that the error estimate for the M^ term is rather conservative. 





AC^) 


input 


P.A. 


z 


UJ 


6 
6 


3/2 
2/3 


-16.83 
-16.84 


-16.73 


7 
7 
7 


4/2 
3/3 
2/4 


-16.89 
-16.91 
-16.90 


-16.49 
-15.85 
-16.79 


8 
8 


4/3 
3/4 


-16.95 
-16.95 


-16.83 



Table 9: Fade results for A^^) 



Mv 





M'i. 


^iw 


M^ln(M^/Af^) 


input 


P.A. 


z 


w 


z 


UJ 


z 


UJ 


6 
6 


[3/2] 
[2/3] 


(*)5.774 
5.340 


— 


— 


— 


-1.21h 
-7.281 


— 


7 
7 
7 


[4/2] 
[3/3] 
[2/4] 


5.197 
5.198 
5.280 


5.338 


11.54 
11.21 


— 


-7.354 
-7.304 
-7.319 


— 



Table 10: Fade results for the coefficients of M^/W^, M^/M^, and M^/Mf 
lnM2/M2,of v4(2). 



Finally we are in the position to write down the decay rate of the top quark up to order 
0{alM^r/Mf). Using the exact results at Born level and at order as in combination with 



Eq. (^ leads to: 



T{t -^ bW) 



a. 



0.8852-2.220—- 15.6(1.1 



TT 



0.788(1) Fo: 



ttc 



TT 



+ ... 



:28) 



where the values Mt = 175 GeV, Mw = 80.4 GeV and as{M^) = 0.11 have been assumed. 
For Mw = the second order QCD corrections amount to roughly 2%. Both at order 
as (see Eq. (p!4D) and a^ the Mw mass corrections "screen" the leading order terms, i.e., 
they arise with negative sign. Whereas the quadratic and quartic corrections at 0{as) 
turn out to be 16% and 3% w.r.t. to the massless result, they amount to roughly 7% and 
1% at 0{al), respectively. 
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Using the 0{al) relation between Mt and the MS mass mt{^) |jTT| one may express the 
result in terms of rrit = mAmt): 



r(t ^ hW) 



M 


mi m^t Tc 


1.28 + 2^^ + 
mf 




( m^\ 
9.16 2 1n ^ * +. 




J\j2 t\jA / 

2.5(8) + 8.1(4) -|: + -f 18.6(5.0)- 

lii't ''H \ 


-4.6(1) In i^+... 


Fo 


0.8576 + 1.98 — + 5.0(1.2) f—) +... 
vr V vr / 






0.9, 


33(1) r 


; 









(29) 

where fo = Gprnc^^Vth^ l{%TX\pX), as{mf) = 0.11, and m^ = 165 GeV has been chosen. At 
0{as) the exact result is used after the second equality sign. As expected, the convergence 
of the perturbative series is better in this case than in Eq. (pS]). 

Let us finally compare the full O(a^) result with the BLM pO[ contributions. They are 
obtained by replacing the number of light fermions in Eq. (|l^) hj —{y~^i) ^^"^ neglecting 
contributions from other colour factors. The result reads: 



^(2) 



3LM 



M^ 



-21.92 + 7.7(8)—^ + 



Mi 



M^ 
Mt 



. s 161 , 
9.9 2 In ^ 

. ^ ' 12 M? 



M! 



(30) 



w , 



For vanishing W boson mass the difference to the complete order a^ result amounts to 
roughly 24% and the (9(M^) term is off by almost 50%. The order of magnitude for the 
quartic term in M^ is reproduced correctly, but the logarithmic term differs from (P7| ) 
by almost a factor of 2. 



4.5 Estimate for b -^ ulu and ji 



ev^Ve 



As it was pointed out in [^, ^, the results for top decay may be used to estimate 
also the QCD corrections to the semi-leptonic decay of the bottom quark. The 0{a1) 
corrections to this process have been obtained recently p2| by computing four-loop on- 
shell diagrams. Nevertheless, using the results of the previous sections we may also derive 
an approximation to this quantity, in this way verifying the consistency of the results of 
2^ , p and the present paper. 

The decay rate for h —>■ ulv can be expressed as 



V{h -^ uW) 



r(o) + 



a. 



TC 



c^^r« + 



rin... 



(31) 



with 



<i) 



2r 



(0) 



dyA^'\y) 



<o) 



Gp 



Vubl^M^ 



192vr3 



(32) 
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where Mf, is the on-shell bottom quark mass and Vub is the CKM matrix element for 6 — >■ m 
transitions. The relation to the top decay rate is established through A^'^\M^r/M^) = A^^\ 
with the A^'^^ defined in Eq. (|13|). Assuming that the functions 

are smooth within < y < 1, one may approximate them by their first few terms in the 
expansion around y = 0: 

^^^Hy) = EK^^+«L?nl0gl/)?/"' (34) 

n>0 

which leads to 

r« = 2rf E(a« /'dyAW(y)y" + ag„ f'dyA^'\y)y-\ogy) . (35) 

^^\ Jo 'Jo J 

For example, in 0*'' approximation, one finds (al^ = for i = 0,1, 2) 

r? = a«r[°) = A«(o)r[°). (36) 

However, A^'^\y) is not really smooth in general. In fact, A'^^\y) has a singularity at 
y = 1 which spoils convergence of the expansion in y. On the other hand, A^'^\y) = 
A^^\y) A^^\y) itself is finite for y = 1. Thus, if a larger number of terms in y is included, 
it is more promising to use directly Eq. (^) and expand the full integrand around small y. 
This is demonstrated in the case of FJ, in Table |lT| where both approaches are compared 
including successively higher powers in y. 

One can see that the approach using Eq. ( P3[ ) provides reasonable estimates for n <, 4, 
where, on the other hand, the results obtained by a naive expansion of the integrand in 
Eq. ( p2[ ) are unsatisfactory. For n > 4, however, the situation becomes opposite: The 
more terms in y are included, the better is the approximation using the latter method. 
The method using Eq. (|35D becomes very unstable. 

The same procedure will now be applied at 0{al). As input we use the numbers of 
Eq. ( p7D and subtract the values for A\ as given in Eqs. (p2D, (pS]), and ([26|), multiplied 



by CpT, according to the transition from n/ = 5 for top decay to n; = 4 for bottom decay. 
Only the central values from these equations will be used, the errors will be suppressed. 



The results are shown in Table |T2|. One can see that both approaches lead to results that 
are fairly consistent with the exact number obtained in P^ . 

Along the same line of reasoning we may derive an estimate for the (9(a^)-corrections 
to the decay rate of the muon. Applying the obvious modifications to the notation of 
Eqs. ( PID and ( P^ and using the results of Eqs. (|2^), (P^D, and (P^D, we find the numbers 
given in Table |13[ Only the method according to Eq. (^) has been applied. Again these 
results agree nicely with the exact results obtained in |]23| . 

This agreement can be considered as a non-trivial check of the results of p2| , p3| and the 
ones obtained in this paper. 
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p(l)/p{0) 


n 


(a) 


(b) 



1 
2 
3 

4 
5 
6 

7 
8 
9 

10 


-2.040 

-1.590 

-2.030 

-1.495 

-2.093 

-1.327 

-2.371 

-0.8998 

-3.071 

0.2045 

-4.881 


-4.080 
-2.580 
-0.3333 
-0.9843 
-1.744 
-1.810 
-1.818 
-1.818 
-1.816 
-1.815 
-1.814 


exact: 


-1.810 



Table 11: Estimates for FJ, using (a) Eq. ([35|), and (b) Eq 



full integrand replaced by its expansion around small y. 
expansion in y that was used as input. 



with the 
n is the order of the 





rl^Vrr 


n 


(a) 


(b) 



1 
2 


-18.6 
-16.8 
-23.4 


-37.1 
-31.2 
-24.9 


exact ||: 


-21.3 



^(2) 



Table 12: Estimates for T[ . Same notation as in Table 11 







n 


77 


elec 


muon 


E 



1 
2 


3.20 
2.38 
4.33 


2.80 
2.50 
3.61 


-0.0636 
-0.0343 
-0.0397 


5.94 
4.85 
7.90 


exact |23(]: 


3.56 


3.22 


-0.0364 


6.74 



Table 13: 



Estimates for r(,^''. 



■'77" denotes the purely photonic corrections, 
"elec" and "muon" the ones involving electron and muon loops, respectively, 
"elec" also includes real emission of an electron-positron pair, of course.) 
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5 Conclusions 

In this work QCD corrections of order a^ to the decay of the top quark into a W boson 
and a bottom quark have been considered. Since the exact treatment of the contributing 
Feynman diagrams is currently out of question, the calculation has been reduced to the 
evaluation of moments. The physical limit is obtained via conformal mapping and Fade 
approximation. The existing results in the limit of a massless W boson could be confirmed 
and new terms of order My^/Mf and M^r/Mf were obtained. Numerically it turns out that 
these power-suppressed terms are rather small. Assuming similar convergence properties 
concerning M-^y for the one- and two-loop corrections we can conclude that the 0{a1) 
corrections for F(t — > Wh) are well under control, including finite PF-mass effects. 

The approach used in this article for the evaluation of the diagrams can certainly be 
carried over to other interesting physical problems, e.g., semileptonic bottom quark decays 
or muon decay as indicated in Section [4.5| . In this paper the reliability of the method has 
been demonstrated by a comparison with a completely different approach. 
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